As discussed in the last meeting, the final model comprises linear
effects for all covariates but Longitude and
Latitude which are included as a two-dimensional additive
effect. This setup leads to the following model output:
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1 ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
## PC2 ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -1.099e+01 1.407e+00 -7.810 5.71e-15 ***
## sand_fraction 2.014e+01 1.092e+00 18.434 < 2e-16 ***
## silt_fraction 1.107e+01 1.850e+00 5.987 2.14e-09 ***
## bulkdensity_soil -2.732e+00 6.977e-01 -3.916 9.00e-05 ***
## ph_soil 1.404e-01 1.753e-01 0.801 0.423134
## soilcarbon -2.469e-01 6.164e-02 -4.006 6.17e-05 ***
## ScenarioSSP1-RCP2.6 -4.811e-01 2.262e-01 -2.127 0.033446 *
## ScenarioSSP3-RCP7.0 1.064e+00 3.323e-01 3.202 0.001365 **
## ScenarioSSP5-RCP8.5 1.508e+00 3.873e-01 3.894 9.85e-05 ***
## time_since_dist 2.642e-04 4.831e-04 0.547 0.584503
## initial_recruitment_BNE 5.121e-03 2.034e-03 2.517 0.011835 *
## initial_recruitment_IBS 2.337e-03 3.400e-03 0.687 0.491794
## initial_recruitment_otherC -1.070e-02 3.924e-03 -2.727 0.006399 **
## initial_recruitment_TeBS 2.879e-02 1.555e-02 1.852 0.064074 .
## initial_recruitment_Tundra -1.725e-03 3.066e-03 -0.562 0.573838
## recruitment_ten_years_BNE -1.013e-04 1.462e-04 -0.693 0.488302
## recruitment_ten_years_IBS 3.509e-04 5.720e-04 0.614 0.539536
## recruitment_ten_years_otherC 2.149e-04 1.184e-03 0.182 0.855931
## recruitment_ten_years_TeBS -4.936e-03 2.733e-03 -1.806 0.070876 .
## recruitment_ten_years_Tundra -1.412e-04 1.981e-04 -0.713 0.475821
## previous_state_BNE 3.901e-02 2.125e-02 1.836 0.066390 .
## previous_state_IBS 4.503e-02 1.568e-02 2.872 0.004073 **
## previous_state_otherC -1.302e-01 3.260e-02 -3.994 6.50e-05 ***
## previous_state_TeBS 1.792e-01 9.542e-02 1.878 0.060389 .
## previous_state_Tundra 1.159e-01 1.284e-01 0.903 0.366787
## PC1_temp 1.112e-01 6.917e-03 16.080 < 2e-16 ***
## PC1_precip -2.051e-04 5.123e-05 -4.003 6.25e-05 ***
## (Intercept).1 4.387e-01 8.405e-01 0.522 0.601736
## sand_fraction.1 5.706e-01 6.517e-01 0.876 0.381271
## silt_fraction.1 -2.358e+00 1.103e+00 -2.138 0.032553 *
## bulkdensity_soil.1 -3.122e-01 4.166e-01 -0.750 0.453542
## ph_soil.1 2.039e-01 1.048e-01 1.945 0.051726 .
## soilcarbon.1 -1.914e-02 3.678e-02 -0.520 0.602827
## ScenarioSSP1-RCP2.6.1 -7.121e-02 1.356e-01 -0.525 0.599382
## ScenarioSSP3-RCP7.0.1 -1.968e-01 1.998e-01 -0.985 0.324606
## ScenarioSSP5-RCP8.5.1 -4.229e-01 2.331e-01 -1.814 0.069697 .
## time_since_dist.1 -9.462e-05 2.877e-04 -0.329 0.742225
## initial_recruitment_BNE.1 1.102e-03 1.212e-03 0.910 0.362992
## initial_recruitment_IBS.1 -7.367e-03 2.025e-03 -3.639 0.000274 ***
## initial_recruitment_otherC.1 -5.035e-03 2.337e-03 -2.154 0.031206 *
## initial_recruitment_TeBS.1 -9.519e-02 9.265e-03 -10.274 < 2e-16 ***
## initial_recruitment_Tundra.1 -4.479e-04 1.826e-03 -0.245 0.806259
## recruitment_ten_years_BNE.1 1.018e-05 8.705e-05 0.117 0.906866
## recruitment_ten_years_IBS.1 1.404e-03 3.407e-04 4.122 3.76e-05 ***
## recruitment_ten_years_otherC.1 9.022e-04 7.051e-04 1.279 0.200740
## recruitment_ten_years_TeBS.1 1.441e-02 1.628e-03 8.851 < 2e-16 ***
## recruitment_ten_years_Tundra.1 -1.334e-04 1.180e-04 -1.131 0.258263
## previous_state_BNE.1 -2.674e-02 1.266e-02 -2.113 0.034598 *
## previous_state_IBS.1 -8.060e-03 9.343e-03 -0.863 0.388337
## previous_state_otherC.1 -8.791e-02 1.942e-02 -4.526 6.01e-06 ***
## previous_state_TeBS.1 -2.610e-01 5.684e-02 -4.592 4.40e-06 ***
## previous_state_Tundra.1 -1.604e-01 7.647e-02 -2.097 0.035966 *
## PC1_temp.1 2.999e-03 4.162e-03 0.721 0.471157
## PC1_precip.1 -1.902e-07 3.079e-05 -0.006 0.995071
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df Chi.sq p-value
## s(Lon,Lat) 26.02 28.40 1023.9 <2e-16 ***
## s.1(Lon,Lat) 26.78 28.66 357.9 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Deviance explained = 66.1%
## -REML = 3280 Scale est. = 1 n = 1443
The effective degree of freedom (edf) of the smooth effects clearly indicates a non-linear relationship, which is confirmed by the visualization of the effects in Figure 1:
Longitude and Latitude for PC 1
(left) and PC 2 (right).As before, looking at the residuals in Figure 2 reveals a structure in the first PC which results from the input data being bounded to the interval \([0,1]\).
This structure also becomes apparent when considering the true PC scores and the predicted ones in Figure 3:
In order to address the structure visible in Figure 3, the input data, i.e., the PC scores need to be transformed. Note that since the second PC does not show any structure, this transformation is only performed for the first PC. The overall shape of the transformation function in Figure 3 resembles to a generalized logistic function:
The formula of this function is given by:
Therefore, the parameters \(A,k,x_0\) and \(C\) need to be set or derived. Figure 5 shows the function with manually derived parameters on the left side. These values were used as starting parameters to determine non-linear least-squares estimates of \(\text{fitted PC 1 scores} \sim \text{glf}(\text{PC 1})\). The resulting curve is shown on the right side of Figure 5.
With that, and the inverse of the generalized logistic function given by
the PC 1 scores are transformed to \(\text{PC1 _trafo} = \text{gls}^{-1}(\text{PC 1})\)
Refitting the previous model with the transformed response yields the following model summary:
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1_trafo ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
## PC2 ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -7.588e+00 1.010e+00 -7.516 5.65e-14 ***
## sand_fraction 1.465e+01 7.836e-01 18.692 < 2e-16 ***
## silt_fraction 7.701e+00 1.327e+00 5.805 6.45e-09 ***
## bulkdensity_soil -1.702e+00 5.006e-01 -3.401 0.000672 ***
## ph_soil -2.135e-02 1.258e-01 -0.170 0.865302
## soilcarbon -1.603e-01 4.422e-02 -3.625 0.000288 ***
## ScenarioSSP1-RCP2.6 -3.928e-01 1.624e-01 -2.418 0.015587 *
## ScenarioSSP3-RCP7.0 7.133e-01 2.387e-01 2.988 0.002808 **
## ScenarioSSP5-RCP8.5 8.881e-01 2.783e-01 3.191 0.001420 **
## time_since_dist 2.192e-04 3.464e-04 0.633 0.526944
## initial_recruitment_BNE 3.203e-03 1.459e-03 2.196 0.028110 *
## initial_recruitment_IBS 6.174e-05 2.438e-03 0.025 0.979795
## initial_recruitment_otherC -7.977e-03 2.814e-03 -2.835 0.004585 **
## initial_recruitment_TeBS 2.172e-02 1.115e-02 1.948 0.051447 .
## initial_recruitment_Tundra 8.040e-04 2.199e-03 0.366 0.714638
## recruitment_ten_years_BNE -8.158e-05 1.048e-04 -0.778 0.436337
## recruitment_ten_years_IBS 1.415e-04 4.102e-04 0.345 0.730181
## recruitment_ten_years_otherC 2.373e-04 8.489e-04 0.280 0.779813
## recruitment_ten_years_TeBS -3.956e-03 1.960e-03 -2.018 0.043540 *
## recruitment_ten_years_Tundra -2.429e-04 1.420e-04 -1.710 0.087323 .
## previous_state_BNE 2.971e-02 1.524e-02 1.949 0.051239 .
## previous_state_IBS 3.671e-02 1.124e-02 3.265 0.001094 **
## previous_state_otherC -6.818e-02 2.338e-02 -2.917 0.003538 **
## previous_state_TeBS 1.319e-01 6.843e-02 1.928 0.053866 .
## previous_state_Tundra 2.079e-01 9.207e-02 2.258 0.023918 *
## PC1_temp 7.947e-02 4.970e-03 15.989 < 2e-16 ***
## PC1_precip -1.796e-04 3.680e-05 -4.880 1.06e-06 ***
## (Intercept).1 4.346e-01 8.405e-01 0.517 0.605122
## sand_fraction.1 5.691e-01 6.517e-01 0.873 0.382560
## silt_fraction.1 -2.359e+00 1.103e+00 -2.139 0.032449 *
## bulkdensity_soil.1 -3.093e-01 4.166e-01 -0.743 0.457753
## ph_soil.1 2.043e-01 1.048e-01 1.949 0.051309 .
## soilcarbon.1 -1.886e-02 3.679e-02 -0.513 0.608228
## ScenarioSSP1-RCP2.6.1 -7.246e-02 1.356e-01 -0.535 0.592986
## ScenarioSSP3-RCP7.0.1 -1.989e-01 1.998e-01 -0.995 0.319573
## ScenarioSSP5-RCP8.5.1 -4.255e-01 2.332e-01 -1.825 0.068006 .
## time_since_dist.1 -9.435e-05 2.877e-04 -0.328 0.742941
## initial_recruitment_BNE.1 1.107e-03 1.212e-03 0.914 0.360729
## initial_recruitment_IBS.1 -7.368e-03 2.025e-03 -3.639 0.000273 ***
## initial_recruitment_otherC.1 -5.040e-03 2.337e-03 -2.157 0.031033 *
## initial_recruitment_TeBS.1 -9.522e-02 9.265e-03 -10.278 < 2e-16 ***
## initial_recruitment_Tundra.1 -4.502e-04 1.826e-03 -0.247 0.805287
## recruitment_ten_years_BNE.1 9.966e-06 8.705e-05 0.114 0.908853
## recruitment_ten_years_IBS.1 1.405e-03 3.407e-04 4.123 3.73e-05 ***
## recruitment_ten_years_otherC.1 9.027e-04 7.051e-04 1.280 0.200479
## recruitment_ten_years_TeBS.1 1.441e-02 1.628e-03 8.854 < 2e-16 ***
## recruitment_ten_years_Tundra.1 -1.332e-04 1.180e-04 -1.129 0.258887
## previous_state_BNE.1 -2.678e-02 1.266e-02 -2.116 0.034314 *
## previous_state_IBS.1 -8.103e-03 9.343e-03 -0.867 0.385824
## previous_state_otherC.1 -8.798e-02 1.942e-02 -4.529 5.92e-06 ***
## previous_state_TeBS.1 -2.612e-01 5.684e-02 -4.595 4.33e-06 ***
## previous_state_Tundra.1 -1.605e-01 7.647e-02 -2.099 0.035858 *
## PC1_temp.1 2.960e-03 4.162e-03 0.711 0.477050
## PC1_precip.1 -6.281e-07 3.079e-05 -0.020 0.983725
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df Chi.sq p-value
## s(Lon,Lat) 26.16 28.46 1183.9 <2e-16 ***
## s.1(Lon,Lat) 26.78 28.66 357.2 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Deviance explained = 67.7%
## -REML = 2815.3 Scale est. = 1 n = 1443
The structure has almost disappeared when looking at the transformed residuals in Figure 6. In particular, it is still present at the limit (\(x=5\)), indicating that the fit of the generalized logistic function is not optimal there.
Looking at the equivalent plot to Figure 5, i.e., true vs. fitted PC scores, confirms that the structure is no longer present.
In order to get an impression of the goodness of fit when considering unseen data, Figure 8 shows the true PC scores plotted against the predicted ones. As in Figure 7, the structure is now no longer present.
In order to get an impression, how the model performs in detail, Figure 9 shows the true functional fit and the predicted one from the model for four randomly picked disturbed grid cells, one for each scenario.
Figure 10 shows the smooth effect of Longitude and
Latitude:
Longitude and Latitude for
PC 1 (left) and PC 2 (right) for the model with transformed PC 1.The interpretation of the model output comprises two steps which need to be kept in mind:
Let’s start with an example. Suppose there is only variable
sand_fraction included in the model and its parameter
estimate is \(\beta_{sand\_fraction} \approx
14.65\), just as for the actual fitted model. The model could
then be rewritten as
with \(\text{glf}()\) being the generalized logistic function. Since this data transformation is non-linear, the interpretation of the parameters is not as straightforward as for a linear model.
Looking at the general structure of the generalized logistic function, one can say that for an increase by one unit it follows that:
This leads to the conclusion that parameter estimates above 0 generally increase the estimated PC 1 scores, while parameter estimates below 0 decrease them.
Considering the whole model output, the derived coefficients can now be classified into different groups for PC 1. Note that significance here means significance on the \(5\%\)-level.
Note that significance here is just based on the model output and does not consider uncertainty in the PC scores.
Figure 11 shows a visualization of the parameter estimates and their confidence intervals:
Since many of the estimates are close to zero, Figure 12 shows the smallest coefficients for a better interpretation:
In order to now draw conclusion on which PFT is mainly influenced by which covariate, Figure 13 shows (again) the first PC:
Let’s go back to the previous example and take a closer look at \(\beta_{sand\_fraction} = 14.65\). Since it is larger than 0, increasing the sand fraction by one unit leads to a higher estimated PC 1 score. Regarding Figure 12, this in turn indicates more Needleleaf evergreen (1), less Pioneering broadleaf (2) as well as more of Conifers (other) (3) and Tundra (5). This PC hardly affects Temperate broadleaf (4) which is covered by the second one.
The interpretation of the other covariates is analogous.
The interpretation of the second PC is less complicated since no transformation of the response is applied. This means in particular, in our example with only one covariate, the model facilitates to
For \(\beta_{sand\_fraction} \approx 0.57\), this means that increasing the sand fraction by one unit increases the predicted PC 2 scores by \(0.57\) on average.
Figure 14 shows a visualization of the estimated coefficients:
Again, Figure 15 shows the smallest effects for a better interpretation:
All coefficient estimates are provided in the following table:
For a valid interpretation of the second PC, Figure 16 shows the second PC for all five PFTS. Again, for variables with estimated parameter values larger than 0, increasing the respective variable by one unit leads to increased PC 2 scores, which in turn indicate a higher share of Needleleaf evergreen (1), Pioneering broadleaf (2) and Conifers (other) (3) in the last six decades of recovery. The second PC especially portrays Temperate broadleaf (4). A unit increase in positive parameter estimates leads to a lower proportion of Temperate broadleaf than on average. For Tundra (5) the effects are only marginal.
Questions here:
So far, four available covariates have not been considered in the
model: tas_yearlymin, tas_yearlymax,
Nuptake and Nuptake_total. All these variables
are functional, and for Nuptake, that is, the nitrogen
uptake per PFT and grid cell, there is even one trajectory for each PFT.
Converting these variables into the required format includes FPCAs for
each of the 5 PFTs. Details on the FPCAs (e.g., plots of the functional
fits, PCs, reconstructions, numbers like accounted variance) are left
out here, but can be presented if needed.
Similar as for tas_yearlymean and
pr_yearly_sum, the first PC resulting of each of the seven
conducted FPCAs are examined in more detail. Note that especially for
the PFT-wise nitrogen uptakes, the first PC accounts for 70-85% of the
variability, so including more PCs might improve the model fit. This is
left out here to not overcomplicate the model.
Figure 17 shows the correlations between all available covariates and the response variable PC 1 and PC 2:
It becomes clear that both PCs are in parts strongly correlated with the PFT-wise nitrogen uptakes. This result is not surprising, since the amount of nitrogen that a PFT is taking is a clear indicator of the share of that PFT.
Question here: It clearly feels like ‘cheating’ using the PFT-wise nitrogen uptake. Should it be removed? Is the correct direction of the causality given here?
The following model output is from a model including only the PFT-wise nitrogen uptakes. The value of deviance explained is already at \(77.5\%\) (!):
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1_trafo ~ PC1_nuptake_BNE + PC1_nuptake_IBS + PC1_nuptake_otherC +
## PC1_nuptake_TeBS + PC1_nuptake_Tundra
## PC2 ~ PC1_nuptake_BNE + PC1_nuptake_IBS + PC1_nuptake_otherC +
## PC1_nuptake_TeBS + PC1_nuptake_Tundra
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.967e-02 3.704e-02 0.531 0.595450
## PC1_nuptake_BNE -4.694e-04 1.045e-03 -0.449 0.653364
## PC1_nuptake_IBS 2.035e-02 3.031e-04 67.138 < 2e-16 ***
## PC1_nuptake_otherC 5.376e-03 1.428e-03 3.764 0.000167 ***
## PC1_nuptake_TeBS 1.615e-03 5.160e-04 3.129 0.001754 **
## PC1_nuptake_Tundra 6.830e-03 9.427e-04 7.245 4.31e-13 ***
## (Intercept).1 7.760e-03 2.487e-02 0.312 0.755010
## PC1_nuptake_BNE.1 -6.080e-04 7.017e-04 -0.867 0.386193
## PC1_nuptake_IBS.1 2.143e-05 2.035e-04 0.105 0.916119
## PC1_nuptake_otherC.1 5.468e-03 9.588e-04 5.703 1.18e-08 ***
## PC1_nuptake_TeBS.1 1.933e-02 3.464e-04 55.798 < 2e-16 ***
## PC1_nuptake_Tundra.1 2.281e-03 6.329e-04 3.605 0.000312 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## Deviance explained = 77.5%
## -REML = 1927.7 Scale est. = 1 n = 1443
Removing the PFT-wise nitrogen uptake and taking the total one instead and all the variables in the initial model with transformed responses also clearly improves the model fit:
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1_trafo ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip + PC1_nuptake_total
## PC2 ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip + PC1_nuptake_total
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -5.026e+00 9.109e-01 -5.518 3.43e-08 ***
## sand_fraction 8.881e+00 7.618e-01 11.657 < 2e-16 ***
## silt_fraction 7.129e+00 1.186e+00 6.013 1.82e-09 ***
## bulkdensity_soil -1.540e+00 4.469e-01 -3.446 0.000568 ***
## ph_soil -2.743e-02 1.122e-01 -0.244 0.806888
## soilcarbon -1.383e-01 3.949e-02 -3.503 0.000461 ***
## ScenarioSSP1-RCP2.6 -3.097e-01 1.448e-01 -2.139 0.032414 *
## ScenarioSSP3-RCP7.0 9.394e-01 2.126e-01 4.419 9.90e-06 ***
## ScenarioSSP5-RCP8.5 1.176e+00 2.478e-01 4.746 2.08e-06 ***
## time_since_dist 2.580e-04 3.097e-04 0.833 0.404737
## initial_recruitment_BNE 3.191e-04 1.313e-03 0.243 0.807994
## initial_recruitment_IBS 3.339e-03 2.186e-03 1.528 0.126606
## initial_recruitment_otherC -5.063e-03 2.520e-03 -2.009 0.044540 *
## initial_recruitment_TeBS 5.716e-02 1.013e-02 5.644 1.66e-08 ***
## initial_recruitment_Tundra 1.775e-03 1.966e-03 0.903 0.366720
## recruitment_ten_years_BNE -6.758e-06 9.378e-05 -0.072 0.942553
## recruitment_ten_years_IBS -4.674e-04 3.681e-04 -1.270 0.204121
## recruitment_ten_years_otherC 4.899e-04 7.589e-04 0.646 0.518535
## recruitment_ten_years_TeBS -9.486e-03 1.775e-03 -5.345 9.03e-08 ***
## recruitment_ten_years_Tundra -3.295e-04 1.271e-04 -2.593 0.009511 **
## previous_state_BNE 3.926e-02 1.363e-02 2.880 0.003977 **
## previous_state_IBS 4.384e-02 1.005e-02 4.360 1.30e-05 ***
## previous_state_otherC -1.877e-02 2.105e-02 -0.892 0.372514
## previous_state_TeBS 2.921e-01 6.172e-02 4.732 2.23e-06 ***
## previous_state_Tundra 3.083e-01 8.246e-02 3.738 0.000185 ***
## PC1_temp 4.466e-02 4.765e-03 9.374 < 2e-16 ***
## PC1_precip -1.148e-04 3.294e-05 -3.487 0.000489 ***
## PC1_nuptake_total 1.354e-02 7.035e-04 19.246 < 2e-16 ***
## (Intercept).1 2.206e+00 7.902e-01 2.791 0.005253 **
## sand_fraction.1 -3.430e+00 6.612e-01 -5.187 2.14e-07 ***
## silt_fraction.1 -2.773e+00 1.029e+00 -2.693 0.007071 **
## bulkdensity_soil.1 -1.680e-01 3.878e-01 -0.433 0.664895
## ph_soil.1 1.856e-01 9.729e-02 1.908 0.056439 .
## soilcarbon.1 -2.901e-03 3.428e-02 -0.085 0.932553
## ScenarioSSP1-RCP2.6.1 2.897e-02 1.254e-01 0.231 0.817290
## ScenarioSSP3-RCP7.0.1 3.254e-02 1.838e-01 0.177 0.859451
## ScenarioSSP5-RCP8.5.1 -1.310e-01 2.141e-01 -0.612 0.540587
## time_since_dist.1 -8.611e-05 2.691e-04 -0.320 0.748979
## initial_recruitment_BNE.1 -9.863e-04 1.141e-03 -0.865 0.387266
## initial_recruitment_IBS.1 -5.168e-03 1.899e-03 -2.721 0.006507 **
## initial_recruitment_otherC.1 -2.862e-03 2.190e-03 -1.307 0.191214
## initial_recruitment_TeBS.1 -7.379e-02 8.798e-03 -8.387 < 2e-16 ***
## initial_recruitment_Tundra.1 4.259e-04 1.709e-03 0.249 0.803164
## recruitment_ten_years_BNE.1 6.939e-05 8.148e-05 0.852 0.394437
## recruitment_ten_years_IBS.1 9.681e-04 3.198e-04 3.027 0.002471 **
## recruitment_ten_years_otherC.1 1.003e-03 6.594e-04 1.521 0.128202
## recruitment_ten_years_TeBS.1 1.092e-02 1.542e-03 7.084 1.40e-12 ***
## recruitment_ten_years_Tundra.1 -2.027e-04 1.104e-04 -1.836 0.066431 .
## previous_state_BNE.1 -1.974e-02 1.185e-02 -1.667 0.095605 .
## previous_state_IBS.1 -1.692e-03 8.733e-03 -0.194 0.846357
## previous_state_otherC.1 -5.493e-02 1.829e-02 -3.004 0.002668 **
## previous_state_TeBS.1 -1.567e-01 5.363e-02 -2.921 0.003484 **
## previous_state_Tundra.1 -9.228e-02 7.165e-02 -1.288 0.197784
## PC1_temp.1 -1.965e-02 4.121e-03 -4.768 1.86e-06 ***
## PC1_precip.1 5.934e-05 2.848e-05 2.083 0.037227 *
## PC1_nuptake_total.1 9.276e-03 6.096e-04 15.217 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df Chi.sq p-value
## s(Lon,Lat) 25.91 28.33 675.6 <2e-16 ***
## s.1(Lon,Lat) 25.45 28.14 505.7 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Deviance explained = 75.8%
## -REML = 2504.4 Scale est. = 1 n = 1443
Now, let’s take a look at the two additional climate covariates, namely the first PCs of the minimum and maximum annual temperature. Including these variables (and excluding all nitrogen related ones) leads to the following model fit:
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1_trafo ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_temp_min + PC1_temp_max + PC1_precip
## PC2 ~ s(Lon, Lat) + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_temp_min + PC1_temp_max + PC1_precip
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -7.173e+00 1.004e+00 -7.143 9.11e-13 ***
## sand_fraction 1.428e+01 7.769e-01 18.382 < 2e-16 ***
## silt_fraction 6.944e+00 1.313e+00 5.288 1.23e-07 ***
## bulkdensity_soil -1.723e+00 4.932e-01 -3.493 0.000477 ***
## ph_soil 1.226e-02 1.241e-01 0.099 0.921310
## soilcarbon -1.552e-01 4.355e-02 -3.565 0.000364 ***
## ScenarioSSP1-RCP2.6 -5.673e-01 1.668e-01 -3.401 0.000672 ***
## ScenarioSSP3-RCP7.0 4.485e-01 2.610e-01 1.718 0.085759 .
## ScenarioSSP5-RCP8.5 6.796e-01 3.004e-01 2.262 0.023682 *
## time_since_dist 2.976e-04 3.431e-04 0.867 0.385771
## initial_recruitment_BNE 2.954e-03 1.444e-03 2.046 0.040772 *
## initial_recruitment_IBS 4.748e-04 2.413e-03 0.197 0.844007
## initial_recruitment_otherC -7.549e-03 2.785e-03 -2.710 0.006723 **
## initial_recruitment_TeBS 2.024e-02 1.104e-02 1.832 0.066906 .
## initial_recruitment_Tundra 1.118e-03 2.176e-03 0.514 0.607326
## recruitment_ten_years_BNE -5.565e-05 1.038e-04 -0.536 0.591879
## recruitment_ten_years_IBS 8.039e-05 4.060e-04 0.198 0.843034
## recruitment_ten_years_otherC 1.404e-04 8.409e-04 0.167 0.867400
## recruitment_ten_years_TeBS -3.736e-03 1.940e-03 -1.926 0.054146 .
## recruitment_ten_years_Tundra -2.787e-04 1.407e-04 -1.982 0.047514 *
## previous_state_BNE 3.061e-02 1.510e-02 2.028 0.042559 *
## previous_state_IBS 3.709e-02 1.112e-02 3.336 0.000851 ***
## previous_state_otherC -6.582e-02 2.313e-02 -2.846 0.004429 **
## previous_state_TeBS 1.279e-01 6.774e-02 1.888 0.059038 .
## previous_state_Tundra 2.128e-01 9.117e-02 2.335 0.019564 *
## PC1_temp 6.005e-02 8.491e-03 7.073 1.52e-12 ***
## PC1_temp_min -1.745e-03 2.792e-03 -0.625 0.531870
## PC1_temp_max 2.974e-02 5.728e-03 5.191 2.09e-07 ***
## PC1_precip -1.229e-04 3.769e-05 -3.262 0.001106 **
## (Intercept).1 5.482e-01 1.667e-01 3.289 0.001005 **
## ScenarioSSP1-RCP2.6.1 -7.116e-02 1.424e-01 -0.500 0.617286
## ScenarioSSP3-RCP7.0.1 -1.662e-01 2.233e-01 -0.744 0.456854
## ScenarioSSP5-RCP8.5.1 -4.101e-01 2.569e-01 -1.597 0.110343
## time_since_dist.1 -1.331e-04 2.905e-04 -0.458 0.646936
## initial_recruitment_BNE.1 1.751e-03 1.217e-03 1.439 0.150217
## initial_recruitment_IBS.1 -8.251e-03 2.035e-03 -4.055 5.01e-05 ***
## initial_recruitment_otherC.1 -4.729e-03 2.360e-03 -2.004 0.045050 *
## initial_recruitment_TeBS.1 -9.572e-02 9.368e-03 -10.218 < 2e-16 ***
## initial_recruitment_Tundra.1 -4.601e-04 1.844e-03 -0.249 0.803014
## recruitment_ten_years_BNE.1 -9.776e-06 8.792e-05 -0.111 0.911462
## recruitment_ten_years_IBS.1 1.583e-03 3.411e-04 4.641 3.47e-06 ***
## recruitment_ten_years_otherC.1 6.403e-04 7.106e-04 0.901 0.367527
## recruitment_ten_years_TeBS.1 1.465e-02 1.644e-03 8.907 < 2e-16 ***
## recruitment_ten_years_Tundra.1 -1.007e-04 1.190e-04 -0.846 0.397382
## previous_state_BNE.1 -2.841e-02 1.279e-02 -2.222 0.026311 *
## previous_state_IBS.1 -1.380e-02 9.284e-03 -1.487 0.137022
## previous_state_otherC.1 -9.393e-02 1.951e-02 -4.816 1.47e-06 ***
## previous_state_TeBS.1 -2.760e-01 5.736e-02 -4.812 1.49e-06 ***
## previous_state_Tundra.1 -2.134e-01 7.625e-02 -2.798 0.005138 **
## PC1_temp.1 9.724e-04 7.333e-03 0.133 0.894513
## PC1_temp_min.1 9.963e-04 2.388e-03 0.417 0.676479
## PC1_temp_max.1 9.373e-05 4.882e-03 0.019 0.984681
## PC1_precip.1 -8.618e-06 3.215e-05 -0.268 0.788648
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df Chi.sq p-value
## s(Lon,Lat) 25.45 28.18 587.5 <2e-16 ***
## s.1(Lon,Lat) 26.62 28.62 300.5 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Deviance explained = 67.9%
## -REML = 2825.5 Scale est. = 1 n = 1443
So including these variables does not lead to substantial differences.
As a final attempt, all new variables are included in the model:
##
## Family: Multivariate normal
## Link function:
##
## Formula:
## PC1_trafo ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
## PC2 ~ s(Lon, Lat) + sand_fraction + silt_fraction + bulkdensity_soil +
## ph_soil + soilcarbon + Scenario + time_since_dist + initial_recruitment_BNE +
## initial_recruitment_IBS + initial_recruitment_otherC + initial_recruitment_TeBS +
## initial_recruitment_Tundra + recruitment_ten_years_BNE +
## recruitment_ten_years_IBS + recruitment_ten_years_otherC +
## recruitment_ten_years_TeBS + recruitment_ten_years_Tundra +
## previous_state_BNE + previous_state_IBS + previous_state_otherC +
## previous_state_TeBS + previous_state_Tundra + PC1_temp +
## PC1_precip
##
## Parametric coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -7.588e+00 1.010e+00 -7.516 5.65e-14 ***
## sand_fraction 1.465e+01 7.836e-01 18.692 < 2e-16 ***
## silt_fraction 7.701e+00 1.327e+00 5.805 6.45e-09 ***
## bulkdensity_soil -1.702e+00 5.006e-01 -3.401 0.000672 ***
## ph_soil -2.135e-02 1.258e-01 -0.170 0.865302
## soilcarbon -1.603e-01 4.422e-02 -3.625 0.000288 ***
## ScenarioSSP1-RCP2.6 -3.928e-01 1.624e-01 -2.418 0.015587 *
## ScenarioSSP3-RCP7.0 7.133e-01 2.387e-01 2.988 0.002808 **
## ScenarioSSP5-RCP8.5 8.881e-01 2.783e-01 3.191 0.001420 **
## time_since_dist 2.192e-04 3.464e-04 0.633 0.526944
## initial_recruitment_BNE 3.203e-03 1.459e-03 2.196 0.028110 *
## initial_recruitment_IBS 6.174e-05 2.438e-03 0.025 0.979795
## initial_recruitment_otherC -7.977e-03 2.814e-03 -2.835 0.004585 **
## initial_recruitment_TeBS 2.172e-02 1.115e-02 1.948 0.051447 .
## initial_recruitment_Tundra 8.040e-04 2.199e-03 0.366 0.714638
## recruitment_ten_years_BNE -8.158e-05 1.048e-04 -0.778 0.436337
## recruitment_ten_years_IBS 1.415e-04 4.102e-04 0.345 0.730181
## recruitment_ten_years_otherC 2.373e-04 8.489e-04 0.280 0.779813
## recruitment_ten_years_TeBS -3.956e-03 1.960e-03 -2.018 0.043540 *
## recruitment_ten_years_Tundra -2.429e-04 1.420e-04 -1.710 0.087323 .
## previous_state_BNE 2.971e-02 1.524e-02 1.949 0.051239 .
## previous_state_IBS 3.671e-02 1.124e-02 3.265 0.001094 **
## previous_state_otherC -6.818e-02 2.338e-02 -2.917 0.003538 **
## previous_state_TeBS 1.319e-01 6.843e-02 1.928 0.053866 .
## previous_state_Tundra 2.079e-01 9.207e-02 2.258 0.023918 *
## PC1_temp 7.947e-02 4.970e-03 15.989 < 2e-16 ***
## PC1_precip -1.796e-04 3.680e-05 -4.880 1.06e-06 ***
## (Intercept).1 4.346e-01 8.405e-01 0.517 0.605122
## sand_fraction.1 5.691e-01 6.517e-01 0.873 0.382560
## silt_fraction.1 -2.359e+00 1.103e+00 -2.139 0.032449 *
## bulkdensity_soil.1 -3.093e-01 4.166e-01 -0.743 0.457753
## ph_soil.1 2.043e-01 1.048e-01 1.949 0.051309 .
## soilcarbon.1 -1.886e-02 3.679e-02 -0.513 0.608228
## ScenarioSSP1-RCP2.6.1 -7.246e-02 1.356e-01 -0.535 0.592986
## ScenarioSSP3-RCP7.0.1 -1.989e-01 1.998e-01 -0.995 0.319573
## ScenarioSSP5-RCP8.5.1 -4.255e-01 2.332e-01 -1.825 0.068006 .
## time_since_dist.1 -9.435e-05 2.877e-04 -0.328 0.742941
## initial_recruitment_BNE.1 1.107e-03 1.212e-03 0.914 0.360729
## initial_recruitment_IBS.1 -7.368e-03 2.025e-03 -3.639 0.000273 ***
## initial_recruitment_otherC.1 -5.040e-03 2.337e-03 -2.157 0.031033 *
## initial_recruitment_TeBS.1 -9.522e-02 9.265e-03 -10.278 < 2e-16 ***
## initial_recruitment_Tundra.1 -4.502e-04 1.826e-03 -0.247 0.805287
## recruitment_ten_years_BNE.1 9.966e-06 8.705e-05 0.114 0.908853
## recruitment_ten_years_IBS.1 1.405e-03 3.407e-04 4.123 3.73e-05 ***
## recruitment_ten_years_otherC.1 9.027e-04 7.051e-04 1.280 0.200479
## recruitment_ten_years_TeBS.1 1.441e-02 1.628e-03 8.854 < 2e-16 ***
## recruitment_ten_years_Tundra.1 -1.332e-04 1.180e-04 -1.129 0.258887
## previous_state_BNE.1 -2.678e-02 1.266e-02 -2.116 0.034314 *
## previous_state_IBS.1 -8.103e-03 9.343e-03 -0.867 0.385824
## previous_state_otherC.1 -8.798e-02 1.942e-02 -4.529 5.92e-06 ***
## previous_state_TeBS.1 -2.612e-01 5.684e-02 -4.595 4.33e-06 ***
## previous_state_Tundra.1 -1.605e-01 7.647e-02 -2.099 0.035858 *
## PC1_temp.1 2.960e-03 4.162e-03 0.711 0.477050
## PC1_precip.1 -6.281e-07 3.079e-05 -0.020 0.983725
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df Chi.sq p-value
## s(Lon,Lat) 26.16 28.46 1183.9 <2e-16 ***
## s.1(Lon,Lat) 26.78 28.66 357.2 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Deviance explained = 67.7%
## -REML = 2815.3 Scale est. = 1 n = 1443
Clearly, according to the deviance explained (\(85.4\%\)) this is the best model, but still, there is the problem with the nitrogen uptake.
Possible next steps: